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ABSTRACT 

We present new predictions for the galaxy three-point correlation function (3PCF) using high- 
resolution dissipationless cosmological simulations of a flat ACDM Universe which resolve galaxy-size 
halos and subhalos. We create realistic mock galaxy catalogs by assigning luminosities and colors 
to dark matter halos and subhalos, and we measure the reduced 3PCF as a function of luminosity 
and color in both real and redshift space. As galaxy luminosity and color are varied, we find small 
differences in the amplitude and shape dependence of the reduced 3PCF, at a level qualitatively con- 
sistent with recent measurements from the SDSS and 2dFGRS. We confirm that discrepancies between 
previous 3PCF measurements can be explained in part by differences in binning choices. We explore 
the degree to which a simple local bias model can fit the simulated 3PCF. The agreement between 
the model predictions and galaxy 3PCF measurements lends further credence to the straightforward 
association of galaxies with CDM halos and subhalos. 

Subject headings: cosmology: large-scale structure of universe — galaxies: formation — galaxies: 
statistics — galaxies: halos 



1. INTRODUCTION 

Observations of the higher-order statistics of the galaxy 
distribution can provide fundamental tests of the stan- 
dard cosmological model. For example, higher-order cor- 
relation functions of the mass are predicted to be zero 
in linear perturbation theory for Gaussian initial condi- 
tions, which are expecte d in the simplest inflation models 
of th e early Univers e ()Peeblesl Il980t iBernardeau et al.1 
l2002t ISzapudil 120051 and references therein). In the 
late Universe, however, non-linear gravitational cluster- 
ing and biased galaxy formation lead to non-Gaussianity 
in the galaxy density field, resulting in non-zero con- 
nected N-point correlation functions (NPCFs) with N > 
2. By studying higher-order galaxy statistics on large 
scales, we can test the nature of the initial conditions; 
on smaller scales, the NPCFs can constrain models of 
biased galaxy formation (e. g., iFrv fc Gaztanagal 119931 : 
iFrieman fc Gaztanagal I1994D and the relationship be- 
tween galaxies and their host dark matter halos. 

The three-point correlation function (3PCF, or (), 
and its Fourier-space equivalent, the bispectrum, are 
the first in the hierarchy of higher-order statistics and 
measure the shape dependence of the number of galaxy 
triplets as a function of scale. The 3PCF is sensitive 
to, for example, the shapes of dark matter halos and 
the presence of filamentary str uctures in the large-scale 
struc ture of the Universe (jSefusatti fc Scoccimarrol 
2005). Since th e pioneering w ork of Peebles and col- 
laborators (s ee iPeeblesI 1 19801 and references therein, 
in particular iGroth fc Peebles! 119771 ). there have been 
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many measurements of the 3PCF and bispectrum 
using a variety of ang ular dPeebles fc GrottJ 1 1975 1: 
Groth fc Peebles! 119771 : IGaztanaga fc Friemanl 1994: 
Frieman fc Gaztanagal 119991: iHuterer et all 1200 It 



Szapudi et all 120011 I2002t iRoss et al l 120061) and red- 
shift (IGaztanaga fc .Friemanl 119941: iJing fc Bornerl 
[Verd~ 



19981: IVerde et al.l 119981: IScoccimarro et al .1 I2001al : 



Fcld man et al.ll2001b Ijing fc Bornerll2004f ) catalogs. 

There has also been considerable theoretical work 
to understand the 3PCF and bispectrum using 
non-linear perturbation theory (see IBernardeau et all 

the 



20021 and refer e nces therein), the halo model (see 
Ma fc Frvl 120001: IScoccimarro et all l2001bl: IWang et all 



20041: iFosalba et al.ll2005l), an d cosmological simu l ations 



dBarriga fc Gaztanagal 120021: IScoccimarro et al.l 119991 : 
IGaztanaga fc Scoccimarroll2005t iHou et al.ll2005f) . 

In recent years, there has been renewed inter- 
est in the 3PCF due to the availability of large 
redshift surveys. These surveys now provide both 
the volume and the number of galaxies required 
to make robust measurements of the 3PCF over a 
range of s c ales. For example , rece nt papers by 
Kavoetal.1 (120041), Ifflkage et all (120051) . I Nichol et all 



20061 ) . iNishimichi et all (|2006[ ). iKulkarni et all (|2007 
have presented measurements of the 3PCF from the 
Sloan Digital Sky Survey (SDSS; lYork et all l2000l) 
as a function of s cale, galaxy luminosity, and color. 
INichol et al.l (|2006l ) also quantified the effect of large- 
scale structures on the shape dep endence of the 3PCF 
Likewise, several recent papers (IJing fc Bornerl 12004 
Baugh et aDl2004l: ICroton et all 120041: IGaztanaga et all 



20051 : iPan fc Sz apudi 20051: ICroton et al.l l2006f ) provide 



new measurements of the 3PCF and high-order correla- 
tions from the Two-degree Fie ld Galaxy Redshift Survey 
(2dFGRS: IColless et al.lL200ll ). 

The main results from these recent SDSS and 2dF- 
GRS analyses of the 3PCF are: i) on large scales, the 
observed 3PCF is in qualitative agreement with expec- 
tations for the growth of structure from Gaussian initial 
conditions; ii) on smaller scales (i.e., in the non-linear 
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and weakly non- linear regimes), the 3PCF measured in 
redshift space scales with the redshift-space two-point 
correlation function (2PCF, or £) as C ~ £ \ consis- 



tent with the "hierarchical clustering" ansatz (|Peeblesl 
1980); Hi) the shape dependence of the reduced 3PCF 
depends at most weakly on galaxy luminosity; iv) the 
amplitude of the 3PCF is larger for elongated triangle 
configurations than for more symmetric triangle shapes 
( Gaztanaga fc Scoccimarrol 120051: i Gaztanaea et al.l l2005; 
iNichol et al.M2006t iKulkarni et al.ll2007ft~ — again consis- 
tent with expectations from non-linear clustering theory 

It is important to make detailed comparisons of these 
3PCF observations with theoretical predictions that in- 
corporate a realistic prescription for modeling galax- 
ies and that account for observational effects such as 
redshift-space distortions. Such comparisons can be car- 
ried out in two ways: either the observations can be 
corrected for the redshift-space effects and compared to 
theory in re al space, e.g., using the projected 3PCF 
(Zhc nsll2004D whic h is analo gous to the projected 2PCF 
(|Zehavi et al.ll20 05), or one can build mock galaxy cata- 
logs from cosmological simulations and measure the the- 
oretical 3PCF directly in redshift space. Since we also 
have access to the real-space 3PCF from such mock cata- 
logs, we can investigate in detail the relationship between 
the real- and redshift-space correlation functions. 

In this paper, we pursue the second of these method- 
ologies, using state-of-the-art high-resolution dissipation- 
less dark matter cosmological simulations. These simu- 
lations have the spatial resolution required to identify 
the dark matter (DM) halos and subhalos that host indi- 
vidual galaxies and at the same time encompass a large 
enough volume to probe large-scale structure in a sta- 
tistically reliable way. The model we use assigns galaxy 
properties (luminosity, color, etc.) directly to these DM 
galactic halos and subhalos using simple, empirically- 
based assumptions about these properties. This ap- 
proach differs in both assumptions and the resolution 
required from methods which build galaxy catalogs by 
statistically assigning several galaxies to each (more mas- 
sive) halo using a Halo Occupat i on Distribut i on ( HOD; 
e.g. iBerlind fc WeinbeTg] 120021 : IWang et all 12004 see 
IKulkarni et alJ (120071 ) for constraints on HOD param- 
eters from the SDSS Luminous Red Galaxy sample's 
3PCF) or from semi-analytic met hods. The method use d 
here was first implemented by iKravtsov et al.l (|2004f) . 
and has been applied successfully to different statist i- 
cal studies, including the 2PCF (IConrov et all | 2006[ ). 
galaxy-galaxy lensing (iTasitsiomi et al.ll2004f ). and close 
pair statistics (|Berrier e t al. 2006 }) among others ( see also 
I Vale fc Ostriker! 12004 12009: IConrov et^[2007h . Here 
we extend the study of this model to the 3PCF as a 
function of luminosity, color, and redshift. Where pos- 
sible, we make direct comparisons in redshift space with 
recent observations. 

In f}2j we describe the simulations used in this paper 
and our methods for constructing mock galaxy catalogs 
based on resolved DM halos. We also review the tech- 
niques used to estimate the NPCFs. In SJ3J we present 
measurements of the 3PCF in both real and redshift 
space for both the dark matter and galaxy catalogs, while 
in 21 we study the dependence of the model 3PCF on 
galaxy luminosity and color. In Sj5j we compare the 
model 3PCF with SDSS observations and discuss the ef- 



fects of binning. We relate the 3PCF to a simple non- 
linear bias model in I n d we summarize and discuss 
our findings. 

2. METHODS 

2.1. Dark matter simulations 

We investigate clustering statistics using cosmological 
iV-body simulations of structure formation in the concor- 
dance, flat ACDM cosmology with ft a = 0.7 = 1 — ft m , 
h = 0.7, and as = 0-9, where fi m , ft a are the present 
matter and vacuum densities in units of the critical den- 
sity, h is the Hubble parameter in units of 100 km s _1 
Mpc -1 , and as specifies the present linear rms mass fluc- 
tuation in spheres of radius 8 /i _1 Mpc. The simulations 
used here were run using t he Adaptive Refineme nt Tree 
TV-body code (ART, see IKravtsov et alJ 119971 for de- 
tails), which implements successive refinements in space 
and time in high-density environments. The primary 
simulation box we use is 120 /i _1 Mpc on a side (hereafter, 
LI 20); the number and mass of each dark matter particle 
are N p = 512 3 « 1.34 x 10 8 and m p = 1.07 x 10 9 h- 1 M Q 
respectively. This simulation has been previously used 
to measure se yeral properties of dark matter halos and 
subhalos (e.g lAllgood et alJ 120061: 1 Wechsler et~all l2006t 
IConrov et all 120061: iBerrier et all 12009) . In order to in- 
clude more massive halos and study the effects of the 
size of the sample on the statistical analysis, we also use 
a second simulation with the same cosmological param- 
eters in a bigger box, with 200 ft, _1 Mpc o n a side (which 
was a lso used to measure halo shapes in lAllgood et al.l 
2006). This box contains N p = 256 3 particles with mass 
m p — 3.98 x 10 10 /i -1 M Q , therefore it will lack low mass 
(and luminosity) objects that are included in the L120 
box. 

From these dark matter samples, virialized concentra- 
tions of particles are identified as halos. In order to find 
these halos and their constituent subhalos (concentra- 
tions of virialized matter inside bigger halos), a variant 
of the Bound Densi ty Maxima halo finding algorithm of 
iKlypin et al.l (1.999) is used. This algorithm assigns den- 
sities to each particle using a smoothing kernel on the 32 
nearest neighbors; centering on the highest-overdensity 
particle, each center is surrounded by a sphere of radius 
ffmd = 50A. kpc. The algorithm removes unbound par- 
ticles when calculating the properties of the halos. The 
halo catalog is complete for halos with more than 50 
particles, which corresponds to a minimum halo mass of 
1.6 x 10 10 hr l M Q for the L120 box and 2.0 x 10 12 h~ l M Q 
for halos in the L200 box. 

Henceforth, we will use the terms "distinct halo" to 
mean any halo that is not within the virial radius of a 
larger halo, "subhalo" to indicate a halo that is within 
the virial radius of a larger halo, and "galactic halo" to 
refer to the halo directly hosting a galaxy. Using this 
terminology, the galactic halo of a satellite galaxy is a 
subhalo while the galactic halo of a central or isolated 
galaxy will be a distinct halo. 

2.2. From halos to 'galaxies ' 

It is expected that galaxy properties depend in detail 
not only on the dark matter clustering, but also on the 
gas dynamics, radiative processes, and feedback mecha- 
nisms that affect the baryonic components. A program 



Modeling galaxy three-point statistics 



to include all those physical processes in building mock 
galaxy catalogs from simulations would require introduc- 
ing a number of free parameters and assumptions which 
could partially obscure the relevant mechanisms that de- 
termine the shape and amplitude of the 3PCF. Our ap- 
proach is instead more empirical — to associate galaxies 
of given properties with simulated dark matter halos and 
subhalos, using halo properties that we can measure in 
the simulation, and to see whether this one-to-one corre- 
spondence predicts a galaxy 3PCF that is consistent with 
the observations. As described below, the assignment of 
galaxies to halos was designed to reproduce certain fea- 
tures of the observed galaxy distribution, but the 3PCF 
was not one of these. As a result, the 3PCF constitutes 
a non-trivial test of this approach. 

The primary galaxy samples used here are created by 
assigning galaxy luminosities and colors drawn from the 
SDSS redshift survey to dark matter halos and subhalos, 
using the maximum circular velocity at z = 0, V max , as 
an indicator of the halo virial mass. V max has been found 
to be a good proxy of the galaxy potential well, which 
is a good indicator of stellar mass (|Kravtsov et al.|[2004t 
iConrov et afll200l . 

A galaxy luminosity in the r— band is assigned to 
each halo by matching the cumulative velocity function 
n(> V max ) of all galactic halos (distinct halos and sub- 
halos) to the observed SDSS r-band luminosity function 
(|Blanton et al.|[2b03b| ) at z — 0.1 (the approximate mean 
redshift of the main spectroscopic SDSS galaxy sample 
used to estimate the luminosity function). To correct 
to z = magnitudes , we use the code kcorrect v3_2 
(|Blanton et al.ll2~003af) . Since the limited size of the box 
gives us an upper limit on the luminosities which can 
be reliably studied (in the statistical sense), and at the 
same time we cannot sample the lowest-luminosity ob- 
jects due to limited spatial resolution, for the L120 box 
we present results for galaxies in the absolute magnitude 
range -19 > M> r l > -22, where M> r l = M r - 51og/i. 
The L200 box has a lower spatial resolution, therefore 
it contains only brighter objects, with M^^S — 20. In 
order to assign co lors to the galax ies, we use the proce- 
dure d escribed bv IWechslerl (|2004[ ) and iTasitsiomi etlrtl 
(|2004t l. This method uses the relation between local 
galaxy density (defined as the distance to the tenth near- 
est neighbor brighter than = —19.7 and within 
cz = 1000 km/s) and color observed in the SDSS, using 
the C MU-Pitt Value Added Catalog constructed from 
DR1 (|Abazaiian et al.ll2003l ). to assign a color to each 
moc k galaxy. We use the d istant observer approxima- 
tion (|Bernardeau et afll2002l ) to obtain positions in red- 
shift space. Table[T]describes the subsamples used in this 
study. 

Although the general association of galaxies with dark 
matter subhalos seems quite robust, the detailed asso- 
ciation of galaxy properties with subhalo properties is 
less clear. While galaxy luminosity is expected to be 
quite tightly connected to velocity, the circular velocity 
of a given subhalo decreases with time due to tidal strip- 
ping as it interacts with its host halo. Observable galaxy 
properties such as luminosity and color will likely be less 
affected by this process. This implies that galaxy ob- 
servables may be more strongly correlated with V maXtacc , 
the maximum circular velocity of the halos at the mo- 
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TABLE 1 

Subsamples at z = 



Box 


Subsample 


Number 
of objects 


Density 
[(fe- 1 Mpc)- 


" 3 ] 


L120 


All objects < -19 


25371 


1.5 x 10" 


2 


L120 


-19 < < -20 


15432 


8.9 x 10" 


3 


L120 


-20 < Nl£ < -21 


8153 


4.7 x 10" 


3 


L120 


red (g — r > 0.7) 


10059 


5.8 x 10- 


■3 


L120 


blue (g - r < 0.7) 


15312 


8.9 x 10- 


■3 


L200 


All objects < -20 


43564 


5.4 x 10- 


3 


L200 


-20 < < -21 


30575 


3.8 x 10- 


■3 


L200 


-21 < < -22 


9921 


1.2 x 10- 


3 


L200 


red (g — r > 0.7) 


17278 


2.2 x 10- 


■3 


L200 


blue (g — r < 0.7) 


23748 


2.9 x 10" 


■3 



ment of accretion onto their host, than with the current 
maximum circular velocity, V maXinow . This conclusion 
is supported by measurements of two-poin t statistics on 
both large and small scal es in simulations (|Conrov et al.l 
I2006t iBerrier et al.ll200(l . Motivated by these consider- 
ations, we construct galaxy catalogs using both V max . acc 
and Vmax,now We also use galaxy catalogs at different 
rcdshifts in order to study the evolution of the 3PCF 
with time. All these additional catalogs use the L120 box 
and have the same spatial density as the sample of halos 
selected by V max . now . We note that the model for as- 
signing color is also not uniquely determined. It is, how- 
ever, sufficient to match the two-point clustering length 
for red and blue ga laxies and several o bserved properties 
of galaxy clusters l)Zehavi et al.ll2005f) . Future measure- 
ments of both clustering statistics and of the properties 
of galaxies in groups and clusters should help to further 
refine these galaxy assignment models. 

2.3. Measuring the 3PCF 

Just as the 2PCF measures the excess probability of 
finding two objects separated by a distance r, the 3PCF 
describes the probability of finding three objects in a 
particular triangle configuration compared to a random 
sample. The probability of finding three objects in three 
arbitrary volumes dVi, dVz, an d dV 3 , at posit ions 7T, r 2 
and r 3 respectively, is given by (|Peebleslll980l ) 

P=[l + £(r 12 ) + £(r 23 ) + £(r 31 ) + ({r 12 , r 23 , r 31 )] x 

n 3 dtW 2 dF 3 , (1) 

where n is the mean density of the objects, ry = 7"j — Tj 
the distance between two objects, £ is the 2PCF, and ( 
is the 3PCF: 

£(r 12 ) = (%i)%2)> (2) 
t(ri2,r23,r 31 ) = (S(r 1 )6(r 2 )S(r 3 )); (3) 

here d is the fractional overdensity in the dark matter 
field or in the distribution of galaxies. Since the 3PCF 
depends on the configuration of the three distances, it 
is sensitive to the 2-D s hapes of the spatial structures , 
at large and small scales (Scfusa tti fc Scoccimarroll2005t 
Gaztahag a fc Scoccimarro ll2005). Motivated by the "hi- 
era rchical" form of t h e N-p oint functions, C 0(1 £, 2 , found 
bv iPeebles fc Grothl (|1975l ). we use the reduced 3PCF 
Q(r, u, a) to present our results: 

Q(r, u,a) — 
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((r,u,a) 



£(ri2)£(r 23 ) + £(r 23 )£(r 3 i) + f (rai)f (r u ) ' 



(4) 



This quantity is useful since Q is found to be close to 
unity over a large range o f scales even th ough £ and C vary 
by orders of magnitude ([Peebles! I1980D . To parametrize 
the triangles for the 3PCF measurements, r = ri2 sets 
the scale size of the triangle, while the shape param- 
eters are given by the ratio of two sides of the trian- 
gle, u = r23/ri2, and the angle between the two sides 
of the triangle a = cos (f 12 • 7^3), where f\2, ?23 
are the unit vectors of the first two sides. Following 
iGaztahaga fc Scoccimarrol (|2005l ). triangles where a is 
close to 0° or 180° are referred to as "elongated config- 
urations" , while those with a ~ 50° — 120° are referred 
to as "rectangular configurations" , 

We calculate the 2PCF using the estimator of 
lLandv fc Szalavl (fl99l . 



£ 



DD - 2DR + RR 
RR ' 



(5) 



Here, DD is the number of data pairs normalized by 
Nd x Njj/2, DR is the number of pairs using data and 
random catalogs normalized by NdNr, and RR is the 
number of random data pairs normaliz ed by Nr x Nr/2. 
The 3 PCF is calculated using the ISzapudi fc Szalavl 
(| 19981 ) estimator: 

DDD -ZDDR+ZDRR- RRR 
C = RRR ' < 6 > 

where DDD, the number of data triplets, is normalized 
by Np/6, and RRR, the random data triplets, is nor- 
malized by N R /6. DDR is normalized by N^Nr/2, and 
DRR by N D N R . 

We estimate the errors using jack-knife re-sampling. 
From each galaxy catalog, we construct sixteen subsam- 
ples of L120 or L200; within each of them we remove a 
different region (30 x 60 2 (/i _1 Mpc) 3 for the L120 box, 
and50xl00 2 (/i^Mpc) 3 for the L200 box) . The variance 
cjk of Q is calculated as: 



°JK 



iV-l " 



N 



i=l 



(Qi - Qif 



(7) 



where N — 16 is the number of subsamples, Qi is the 
value for the i— th subsample, and Qi is the mean of the 
Qi. We note in passing that the validity of jack-knife 
resampling as a method to estimate the errors has not 
been explicitly tested with mock catalogs for three-point 
statistics. Although it is beyond the scope of this paper, 
this may be an interesting topic of future investigation 
especially once the statistical power of the measurements 
improves. 

To compute the 2PCF and 3PCF, we use the NPT 
software developed in collaboration with the Auton Lab 
at Carnegie Mellon University. NPT is a fast implemen- 
tation of the NPCFs using multi-resolution kd-trees to 
compute the number of pairs and triplets in a dataset. 
For more de t ails a n d information on the a lgorithm, see 
iMoore et all (|2001| ). lGrav et alJ (|2004D . and lNichol et all 
(2001). 
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We estimate the reduced 3PCF for the distribution of 
dark matter and for galaxies for different triangle con- 
figurations, focusing on the scale and shape dependence 
of the 3PCF. We also investigate its time evolution and 
how it depends on the selection criterion for subhalos. 
We study Q(r, it, a) in both real and redshift space, in 
order to compare our results with current observations 
and disentangle galaxy biasing effects from those which 
are consequences of redshift distortions. 

In order to distinguish scale and, most importantly, 
shape effects, and to keep the errors as small as possi- 
ble, we have chosen an intermediate-resolution binning 
scheme. For studies of equilateral triangles (u = 1 and 
a = 7r/3 rad), we use bins of size Alog(r) = 0.1. For 
measurements of the shape dependence of the 3PCF, we 
use triangles with four different scales r = 0.75, 1.5, 3, 6, 
and 9 /i _1 Mpc, using the shape parameters u = 2, and 15 
angular bins separated by Aa = tt/15 rad; the resolution 
of the bins is given by Ar*jj = ±0.03r,j . This resolution is 
sufficient to see the most important features of the 3PCF 
even on small scales, although it is not sufficient to dis- 
tinguish the "finger-of-god" effect at the smallest scales 
in redshift space, where Q(a) varies very little except 
at very small or elongated a ngles, where it increases to 
many times the mean value (jGaztanaga fc Scoccimarrol 

3.1. The 3PCF in real space 

Figures Q] and [2] show the 3PCF for dark matter par- 
ticles and for galaxies from the L120 and L200 simula- 
tions. Here we plot results for dark matter (thick solid 
line), galaxies in halos selected by V ma x,now (thin solid 
line), and galaxies in halos selected by V maXyacc (long- 
dashed line) for z = and for halos selected by V max ,now 
at z = 1 (short-dashed line) and z = 2 (dotted line). 

Figure ft] shows the reduced 3PCF for equilateral trian- 
gles, Q eq (r), in real (left panel) and redshift space (right 
panel). In real space, the reduced 3PCF for the dark 
matter is only weakly scale dependent on small scales, 
decreases rapidly with increasing scale around r ~ 3/i _1 
Mpc, and falls off more slowly on larger scales. This 
behavior is broadly consistent wit h previous N-body re- 
sults (e.g. lScoccimarro et al.lll998t ) and with expectations 
from leading order non-linear perturbation theory on the 
largest scales (shown as the thick red long-dashed curve 
in Figure 1 left panel), with loop-corrected perturbation 
theory on intermediate scales where the rms perturba- 
tion amplitude S(r) is of order unity (the transition to 
the strongly non-linear regime; adding more orders to the 
calculation would increase the agreement to the iV-body 
dark matter 3PCF amplitude), and with quasi-stable hi- 
erarchical clustering on the smallest scales. Tests with 
the L200 box indicate that the downturn in Q eq for dark 
matter at scales larger than r ~ 8/i _1 Mpc is likely due 
to finite volume effects. 

At scales below r ~ lO/i -1 Mpc, the dark matter Q eq is 
larger than that for the galaxies; this behavior is broadly 
expected if galaxies are more strongly clustered than 
(positively biased with respect to) the mass, cf. eqn. (fTO]) . 
At higher redshift, evolution is seen in Q eg (r) that is con- 
sistent with expectations from non-linear gravitational 
evolution: on the largest scales, the amplitude of Q eq is 
unchanging, as predicted from leading order perturbation 
theory, while the sharp break associated with the tran- 
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real space, L120 box 

dark matter, z-0 

halos, z-0 

- halos, z=l 
halos, z=2 
halos, z=0. V m 
PT, dark matter 




redshift space, L120 box 




r (Mpc/h) 



r (Mpc/h) 



Fig. 1. — The reduced 3PCF, Qeq(r), for equilateral triangles in the L120 box. Left: Results in real space. Thick solid line (black): dark 



matter; thin solid line (black): galactic halos selected by Vm 



at z = 0; short-dashed (green): galactic halos selected by V m ax,now at 



1; dotted (red): galactic halos selected by Vmax,now at z = 2; long-dashed (blue): galactic halos selected by V m ax,acc at z = 0; thick 
long-dashed (red): leading-order perturbation theory, dark matter. Right: Results in redshift space. Line types correspond to the same 
dark matter and halo samples as in the left-hand plot. Error bars are calculated using jack-knife resampling and are only shown for one of 
the samples for clarity. 



sition to the strongly non-linear regime moves to larger 
scales as the density perturbation amplitude increases 
with time. 

Comparing results for subhalos selected by V max ^ now 
and by V max ^ acc , differences in the amplitude of Q eq ap- 
pear on small scales, r < 3 hr 1 Mpc. In halo-model 
language, on these scales the 3PCF is sensitive to the 
internal structure of halos, i.e., to the one- and two-halo 
terms, while the three-halo term dominates the 3PCF on 
larger scales (|Wang et all 120041: iTakada fc Jainll2003D . 

In Figure^ the top panels show how the reduced 3PCF 
depends on triangle shape in real space. In general, 
the 3PCF for elongated configurations is greater than 
for rectangular configurations. This is a consequence of 
the fact that, in non-linear gravitational instability, ve- 
locity fl£WB_te2d_to_oroijT_along gradients of the density 
field (jBernardeau et al J 1200^ 1. The 3PCF is larger for 
dark matter than for galaxies for all shapes and scales, 
although the difference is larger for elongated configura- 
tions. The difference in 3PCF amplitude between rect- 
angular and elongated configurations is larger on large 
scales, in broa d agreement with leadin g-order theoreti- 
cal predictions (jBernardeau et al]l2002f ): on large scales, 
the strong shape dependence is determined by pertur- 
bative non-linear dynamics; on smaller scales, the shape 
dependence is washed out since the coherence between 
the velocity and density fields gives way to virialized mo- 
tions. This scale dependence of the 3PCF shape is also 
reflected in the redshift evolution: in the L120 box (left 
panel), the galaxy 3PCFs at z = 1 and 2 (green-dashed 
and red-dotted curves) essentially retain the primordial 
shape dependence of leading-order non-linear perturba- 
tion theory, i.e., at those redshifts, these scales are still 
close to the quasi-linear regime. At r = 3 h~ 1 Mpc, the 
largest evolution in Q(a) is found for elongated configu- 
rations. 

As was seen in Figure [TJ the effect of changing halo se- 
lection from V ma x,acc to V ma x.now on the 3PCF shape ap- 



pears only on small scales, r < 1.5/i -1 Mpc, i.e., roughly 
within the scale of a typical cluster-mass host halo. 

On the larger scales probed in the L200 box (right 
panel of Figure [2]), the galaxy reduced 3PCF (open cir- 
cles) tracks the shape of the dark matter 3PCF fairly 
well. The difference between the galaxy and dark matter 
3PCF amplitudes on these scales is reasonably well fit 
by a simple bias prescription: the thin blue curve is the 
biased 3PCF that results from fitting the galaxy 3PCF 
with eqn. (TO}; see We also see that the jack-knife 
errors increase on the largest scales, where the effects of 
the finite box size start to become evident. For compar- 
ison, the red long-dashed curve is the 3PCF of the dark 
matt er from leading-order non-linear perturbatio n the- 
ory |li7HMdi^iLan Ull U^A3m^^J^- On 
the largest scales, it is in reasonable agreement with the 
measured 3PCF for the dark matter. 

3.2. The 3PCF in redshift space 

Redshift distortions have been studied in depth (and 
are a useful tool to constrain c osmological parameters) 
for the power spectrum (e.g., iBernardeau et al.1 120021 : 
Ida Angela et all 120051: iTinkeri l2007ft and for th e bis- 
pectrum (IScoccimarro et al. □ 19991 IVerde et all l200l 
iSefusatti et alj I2006T) . Some comp arisons have been 
made for the 3PCF as w ell dMatsubara fc Sutol 
1991 ITakada fc Jainl J200l IWang et al.l 12004ft . 
Gaztanag a fc Scoccimarrol (|2005h found that the 



redshift distortions do not have a strong dependence on 
the cosmological parameters. 

The right panel of Figure [TJ and the bottom panels 
of Figure [2] show the 3PCF in redshift space. The first 
feature that can be seen is a dramatic decrease in the 
amplitude and in the scale and shape dependence of Q 
compared to the real-space measurements. For exam- 
ple, for equilateral triangles, the redshift space Q z (r) 
is reduced compared to the real space Q(r) at small 
scales and increased with respect to the real-space re- 
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Fig. 2. — Measurement of the reduced 3PCF as a function of triangle shape, Q(a), for different scales r, with side ratio u = r2s/ri2 = 2 
fixed. Left: Results in the L120 box for r = 0.75,1.5, and 3 h~ x Mpc (from left to right) in real space (top panels) and redshift space 
(bottom panels) for dark matter and galaxies; line types correspond to the same dark matter and halo samples as in Figure 1. Right: 
Results in the L200 box for r = 6 (left) and 9 h — 1 Mpc (right) in real (top) and redshift space (bottom); Thick solid line (black): dark 
matter; open circles: galaxies in halos selected by V m ax,now at z = 0; thin solid line (blue): predicted galaxy 3PCF using the dark matter 
3PCF 

Qdm and eqn. H10I I with best-fit bias parameters ci,C2 obtained from fitting the galaxy 3PCF at r — 9h ^ Mpc (see Sj6l ; long 
dashed (red): leading order non-linear perturbation theory prediction for dark matter reduced 3PCF. Error bars calculated using jack-knife 
resampling method. 

suits on larger scales. The overall effect is that Q z {r) 
is nearly scale-independent, i.e., the clustering appears 
more hierarchical in redshift space (|Suto fc Matsubar 
11994 iMatsubara k. Sutolll994t IScoccimarro . et al.lll99 ' 
Moreover, in redshift space, the suppression of the galaxy 
3PCF relative to that of the dark matter is apparent 
on all scales; it appears to be relatively independent of 
scale and configuration and is larger than the relative 
suppression in real space, consistent with earlier results 
(|Wang et all 12004 iGaztanaga et al.ll2005f ). It also ap- 
pears that there is very little redshift evolution of the 
galaxy 3PCF in redshift space; the measurements at 
2 = 1 and z = 2 are nearly indistinguishable from each 
other. With regard to halo selection, as in real space we 
find that Qvmax.now < Qvmax.acc, but the differences 
between them are smaller than in real space. 

Together, these results suggest that the shape and scale 
dependence of the reduced 3PCF in redshift space on the 
scales shown here are largely determined by redshift dis- 
tortions, with non- linear gravitational evolution playing 
a subdominant role. 
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4. OBSERVING THE 3PCF: LUMINOSITY AND COLOR 
DEPENDENCE 

To investigate the dependence of the three-point clus- 
tering on galaxy luminosity and color and to make direct 
comparisons with measurements from recent redshift sur- 
veys, we calculate the 3PCF for galaxies with luminos- 
ity and color cuts similar to those that have been ap- 
plied to redshift survey data samples, with luminosity 
and color information obtained as described in §2. The 
luminosities are assigned according to V m ax,now m the 
L120 and L200 boxes, since we hace those measurements 
for both boxes. As seen in the previous section, the use 
of V m ax,now is justified since the differences with respect 
to V max ,acc are only significant on small scales (r^S 1.5 



Fig. 3. — The reduced 3PCF, Q(r), for equilateral triangles as a 
function of galaxy luminosity and color in real and redshift space. 
Top left: Q eq (r) in real space for galaxies divided into luminosity 
bins; long dash-dotted (cyan): —19 > > —20; short-dashed 

(green): -20 > Rl£ > -21; dotted (magenta): -21 > > -22. 
The brightest sample comes from the L200 box, the other two are 
from the L120 box. Bottom left: Results in redshift space; line 
types correspond to the same galaxy samples as in the top left 
panel. Top right: Q zq (r) in real space for galaxies divided accord- 
ing to color, using the L120 box; long-dashed (red): red galaxies 
(g — r > 0.7); short dash-dotted (blue): blue galaxies (g — r < 0.7). 
Bottom right Results in redshift space; line types correspond to 
the same galaxy samples as in the top right panel. Error bars are 
calculated using jack-knife respampling. 

/i -1 Mpc) in real space and are almost completely dis- 
sapear on redshift space. We adopt the same binning 
scheme used in the previous section. 
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Fig. 4. — Q(ct) shape dependence for galaxies divided by lu- 
minosity, for r = 1.5, 3, 6, and 9 h _1 Mpc, with fixed u = 2 in 
real (lines) and redshift space (symbols). Top two plots are re- 
sults for the L120 box: filled squares and long dash-dotted (cyan): 
— 19 > M£ > —20; filled triangles and short-dashed (green): 
-20 > > —21; triangles are slightly shifted to the right for 
clarity. Bottom plots show results for the L200 box: filled trian- 
gles and short-dashed (green): —20 > M™ > —21; filled circles and 
dotted (magenta): —21 > > —22; circles are slightly shifted 
to the right for clarity. Error bars are calculated using jack-knife 
resampling and are shown only for one of the samples. 

4.1. Luminosity dependence 

The left panels of Figure [3] show results for the re- 
duced 3PCF for equilateral configuratons in two lumi- 
nosity bins, in real and redshift space. For equilateral 
triangles, there is a small difference in 3PCF between 
the luminosity samples in real space: the fainter galaxies 
have larger Q(r) than the brighter ones, as expected in 
a simple linear bias model (eqn. [T0|) . The redshift-space 
3PCFs for these galaxies are roughly constant with r, 
Qz,gal( r ) ~ 0.7, in agreement with SDSS m easurements 
for th ese configurations (see Figures 7-9 in iKavo et al.l 
•2(101;. There is a very slight difference between the re- 
duced 3PCF amplitudes for different luminosity bins in 
redshift space, but it is not statistically significant for a 
dataset of this size. This result q ualitatively agrees with 
the SDSS results of IKavo et alj (|2004j ) who also found 
almost no luminosity dependence of the reduced 3PCF 
on these scales. They found a slightly higher amplitude 
for the reduced 3PCF for the -19 > > -20 sample 
compared to that of brighter galaxies, but their results 
for the two luminosity samples were consistent within 
the error bars. The brightest galaxies (dotted line) show 
significant fluctuations in both real and redshift space; 
we think this behavior is due to the small size of the box 
and the low density of objects. 

Figure 2] shows the dependence of Q(a) on galaxy lu- 
minosity in real (curves) and redshift space (points). We 
use the same ordinate scale for all the plots to empha- 
size where configuration effects are more important. The 
top panels show results for small scales, calculated with 
the L120 box. On these scales, the reduced redshift- 
space 3PCF of the brighter sample (-20 > > -21, 
short dashed curve for real space, filled triangles for red- 
shift space) is slightly lower than for the fainter sample 



Fig. 5. — Q(c) shape dependence for galaxies divided by color, 
in the L120 box (top) and L200 box (bottom) in real (lines) and 
redshift space (symbols). Open squares and long- dashed line (red): 
red (g — r > 0.7) galaxies, open triangles and short dash-dotted 
(blue): blue (g — r < 0.7) galaxies; squares are slightly shifted 
to the right for clarity. Error bars are calculated using jack-knife 
resampling and are shown only for one of the samples. 

(-19 > > -20, long dashed-dotted curve for real 
space, filled squares for redshift space) for all angles, 
consistent with the results for Q z (r) using equilateral 
triangles. The lower plots in Figure 2] show the lumi- 
nosity dependence on larger scales, measured using the 
L200 box; note that the luminosity bins in the lower 
plots are -20 > > -21 and -21 > Mz[ l > -22. 
The characteristic U-shape of Q(a) appears clearly in 
the redshift-space measurements on scales larger than 
r = 6 /i _1 Mpc. The luminosity dependence of Q on 
these scales appears non-existent in redshift space and 
only slight in real space. There are hints that the reduced 
3PCF for fainter galaxies may have slightly higher am- 
plitude and stronger shape dependence than for brighter 
galaxies, but these trends are not statistically significant 
in the samples studied here. Rather, the strong lumi- 
nosity dependence observed for the 2PCF (jZehavi et al.l 
2005) appears to be closely matched by a correponding 
dependence of the 3PCF, such that the reduced 3PCF Q 
is roughly independent of luminosity. 

4.2. Color dependence 

The right panels of Figure [3] show the reduced 3PCF 
Q{r) for equilateral triangles separately for red and blue 
galaxies in real (top) and redshift space (bottom). Due 
to limited statistics, we use the full lum i nosity range in 
each color bin, as opposed to lKavo et al.l (|2004D . who di- 
vided the color bins into luminosity subsamples as well 
(see their Figure 9). Since the luminosity functions for 
red and blue galaxies differ, our red and blue samples 
have different characteristic luminosities — the red sam- 
ple is on average brighter. In both real and redshift space, 
the red galaxies appear to have a slightly higher reduced 
3PCF amplitude. Comparing with the left panels, this 
difference is in the opposite sense from that expected 
from the fact that the red galaxies are brighter; put an- 
other way, if we were to compare red and blue samples 
of the same luminosity, the color difference in Q would 
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likely be larger than that seen here. On the other hand, 
we should not overinterpret these tren ds, since the differ- 
ences are within the statistical errors. iKavo et all (120041 ) 
found a similarly weak dependence of Q in redshift space 
on color. 

Figure [5] shows the color dependence of Q(a) in real 
(curves) and redshift space (points), for the same sam- 
ples and configurations as in Figure [4] In real space, on 
scales larger than about 3 hr 1 Mpc, Q for red galaxies 
is larger for rectangular configurations and smaller for 
elongated configurations than for blue galaxies. This be- 
havior is qualitatively consistent with a picture in which 
red galaxies preferentially occupy the inner regions of 
clusters, while blue galaxies tend to trace out more el- 
liptical or filamentary structures. In redshift space, the 
trend with color is largely washed out. These redshift 
space results appear more consistent with t he observa- 
tions of the 2dFGRS (|Gaztanaga et al.ll2005f) . where the 
color differences for Q in r edshift space are s maller than 
those seen in the SDSS bv IKavo et all (l200l . 

The fact that red galaxies have a larger two-point clus- 
tering amplitude than blue galaxies and that the reduced 
3PCF for red galaxies is also larger than for blue galax- 
ies (at least for rectangular configurations) suggests that 
red galaxies have a larger quadratic (non-linear) bias, 
Cf. eqn. [lOl This is qualitatively consistent with the ob- 
served morphology-density or color-density relation ac- 
cording to which red galaxies are preferentially found in 
dense regions, since the latter contribute more strongly 
to the higher-order correlations. The qualitative agree- 
ment between the two- and three-point observations and 
our method for assigning galaxy colors confirms that the 
color of a galaxy, a consequence of many physical pro- 
cesses occurring inside the galaxies, depends largely on 
the surrounding environment. However, we note that 
the color assignment for these subhalos is the most un- 
certain part of the model; future work will be required 
to determine how clustering statistics will change if more 
sophisticated schemes are adopted. 

4.3. Redshift distortions: Galaxy Type and Evolution 

So far, we have explored some differences between mea- 
surements in real and redshift space of the 3PCF for dif- 
ferent simulated galaxy samples. Here we investigate in 
more detail whether the redshift distortions of the 3PCF 
are universal or instead depend on the type of galaxy 
studied and to what extent they evolve with time. We 
study the behavior of the quantity 

AQ(r, u, a) = Q z (r, u, a) - Q r (r, u, a), (8) 

where Q r , Q z represent the reduced 3PCF measured in 
real and redshift space, respectively. We explore AQ 
as a function of galaxy type (luminosity and color) and 
epoch. Figure [6] shows AQ for equilateral triangles as a 
function of scale and also the configuration dependence 
for triangles with fixed r, u = 2, and different opening 
angles a. 

In general, the trends are similar to those seen above: 
Qz < Qr for small scales and for elongated triangle con- 
figurations, while the opposite behavior is seen for larger 
scales and for rectangular configurations. At z = 0, for 
equilateral triangles AQ(r) appears to display a roughly 
universal scale dependence, independent of galaxy type, 
with AQ(r) ~ 0.67r~ 2 - 2.6r~ 1 - 0.02r + 0.88 over the 
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Fig. 6. — AQ(r, u, a) = Qz (r, u, a) — Qr(r, u, a) as a function of 
galaxy type and redshift. Top: AQ for dark matter at z = and for 
galaxy samples at different redshifts in the L120 box, for equilateral 
triangles (top left) and as a function of angle for triangles with 
r = 3 h~ 1 Mpc and u = 2 (top right). Bottom: AQ as a function 
of galaxy type for equilateral triangles in the L120 box (bottom 
left) and as a function of angle for r = 6 h _1 Mpc and u = 2 in the 
L200 box (bottom right). 

range 0.5 < r < 5 /i _1 Mpc. On the other hand, for 
r = 6/i _1 Mpc and u — 2, the shape dependence of 
AQ shows more dependence on galaxy type, with blue 
galaxies having larger values than red galaxies and bright 
galaxies smaller values than faint galaxies. Note that the 
redshift distortions of Q appear insensitive to whether 
the subhalos are identified at the present or at the time 
they are first accreted onto host halos. 

We see clear evolution of AQ with redshift, as expected 
since redshift distortion effects should become more pro- 
nounced as perturbations become more non-linear. For 
equilateral triangles, the upper left panel of Figure [6] 
shows that the scale r where AQ ~ increases with 
time. The shape dependence of AQ at fixed scale also 
appears to increase with time. 

5. OBSERVING THE 3PCF: COMPARISON TO SDSS DATA 
AND BINNING EFFECTS 

Below, we compare our res ults with mea s urem ents of 
the 3PCF from the SDSS bv iNichol et~ai] (pOOl . The 
SDSS sample is magnitude limited, with m r < 17.5, and 
has additional cuts in absolute magnitude, — 19 < M r < 
—22, and in redshift, 0.05 < z < 0.15. In order to com- 
pare with this data, we randomly resample the galaxies 
from the L120 simulation (in particular, the V max ^ noew 
sample at z — 0), so that they have the same absolute 
magnit ude distribution as t he SDSS flux-limited sample 
used in INichol et alJ (|2006i ) ; we call these pseudo-flux- 
limited samples. This resampling technique does not 
properly model the distance-dependent selection func- 
tion of the SDSS, but it should reproduce its clustering 
properties on average. Moreover, since, as we have just 
shown, the luminosity dependence of the reduced 3PCF 
is weak, we expect that this procedure should be suffi- 
ci ently accurate for o ur purposes. 

INichol et all (2006) measure the shape dependence of 
Q(a) for four different length scales, using wide bins in 
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Fig. 7.— Top: Q z (r,u = 2, a) at r = 1.0, 2.0 and 4.0 fa" 1 Mpc 
from SDSS observations (green points) from[Nichol ct al. (120061 ) 
and from our pseudo-fiux-limited-sample in the L120 box (long- 
dashed lines) with the same binning scheme used in the mentioned 
paper. Bottom: Effect of binning in the 3PCF measurements. 
Q z (r, u = 2, a) at r = 2.0 and 4.0 h _1 Mpc in the L120 box and 
r = 10 /i — *Mpc in the L200 box for the simulations in a volumc- 
limited sample, so lid (red): results using a wide binning scheme 
UNichol ct al. 2006); long-dashed (black): results with narrow bin- 
ning. In the middle panel, we plot the results using SDSS data 
with wide (green points) and narrow (triangles) binning. 



r,u, and a (see their paper for more details): Ar = 1 
/i _1 Mpc, Au = 1 and Aa = 0.1 rad. In Figure El 
we compare the SDSS results (green points) to the cal- 
culation of Q z (r,u — 2, a) using our model redshift- 
space, volume-limited sample 3PCF (solid lines) and 
our pseudo-magnitude- limited sample (dotted line), each 
with the same binning as the data. We use the L120 box 
for small scales and the L200 box for the r = 10 /i _1 Mpc 
measurements. As the top panels of Figure [7] show, the 
model and data show good agreement within the model 
jack-knife error bars in amplitude as well as in the shape 
ofQ(a). 

Over the range of scales shown here, both the SDSS 
data and the simulations show relatively little shape de- 
pendence for the reduced 3PCF amplitude. This is in 
contras t to the results of j pfl abo ve and to the 2dFGRS re- 
sults of IGaztafiaga et all ( 20051 ) on similar scales, where 
a significantly stro nger shape dependence of Q(a) is evi- 
dent. As noted by IGaztafiaga fc Scoccimarrol (|2005l ) and 
iKulkarni et alj (|2007f ). the differences can be traced to 
the binning scheme: the relatively wide binning scheme 
used here results in smearing and therefore suppression 
of the U-shape of Q over most scales of interest. This 
effect is not due primarily to the binning in a: small 
bins in both r and u are necessary to see the effects of 
shape-dependent clustering. 

To illustrate these effects in more detail, in the lower 
panels of Figure [7] we also show results for the same sam- 
ple but with a narrower binning scheme, using Ar = 
0.1 /i _1 Mpc (ten times smaller), Au = 0.2 (five times 
smaller), and Aa = 0.05 rad (two times smaller). With 
the narrower bins (shown by the dashed curves in Figure 
[7|) , the shape-dependence of Q is more pronounced for the 
simulation, especially on larger scales. For comparison, 



Fig. 8. — Galaxy vs. dark matter overdensities for all galaxies 
in the L200 box measured in randomly placed spheres of radius 
XQh~ 1 Mpc in real (left panel) and redshift space (right panel). 
Solid black curve denotes best fit using the quadratic bias relation 
of eqn. [9] Long-dashed red line indicates the best linear bias fit 
(&2 = 0) to the 2PCF, i.e., estimati ng = (, ga l/(,dm,- Short-dashed 
blue curve indicates best fit of eqn. riOTto the 3PCF of the galaxies 
and dark matter. 

for r — 4 /i _1 Mpc (lower center panel), the solid trian- 
gles shows results for the SDSS flu x-limited sample us - 
ing a similar narrow binning scheme iNichol et al. (2006), 
again showing good agreement between the model and 
the data. 

6. GALAXY BIAS AND THE 3PCF 

As we have seen, the 3PCF predicted for galax- 
ies differs systematically from that expected for dark 
matter. These differences reflect differences in the 
spatial distributions of these two populations; higher- 
order statistics can therefore provide important con- 
stra ints upon the bias between galaxies and dark mat- 
ter (|Frv fc Gaztanagdl993tlFrieman fc Gaztanagall994f ) 
and its dependence upon galaxy properties. 

On large scales, where the rms dark matter and galaxy 
overdensities are small compared to unity it is com- 
mon to adopt a deter ministic, local bias model (e.g. 
iFrv fc Gaztanagalll993l ). 



dvn ) 



hS. 



dm 



2 



(9) 



where S ga i and 5dm are the local galaxy and dark matter 
overdensities smoothed over some scale R. We can use 
the simulations above to test how well this simple bias 
prescription characterizes the galaxy distribution and its 
clustering statistics. 

Figure [8] shows the relation between S ga i and 5dm for 
all subhalos in the L200 box. The points show the over- 
densities for the galaxy and dark matter fields in ran- 
domly placed spheres of radius 10/i _1 Mpc in both real 
(left panel) and redshift space (right panel). The solid 
black curves show the best quadratic fits of the form 
in eqn. @. The quadratic local bias model appears 
to do a reasonable job in characterizing the mean rela- 
tion. Nonetheless, there is significant scatter, either due 
to stochastic bias or dependence of bias on other prop- 
erties than 5dm, that is not captured by this simple bias 
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model. The errors on these fits, also extended to samples 
divided by galaxy luminosity, are shown in Fig. [9] by the 
light solid contours. 

To leading order, this bias prescription leads to a 
relation between the galaxy a nd dark matter reduced 
3PCF amplitudes of the form (iFrv fe Gaztanagal ll993t 
iGaztanaga fc Scoccimarro1l2005f ) 

Qgal = — (Qdm + Cfc) , (10) 
^1 

where ci = b\ and c 2 = bi/bx. Also to leading order, 
at low overdensities the relation between the galaxy and 
dark matter 2PCF amplitudes in this model is given just 
by the linear bias, £ ga i = b\^dm- We can test how well 
this bias prescription captures the clustering statistics 
by fitting these relations to the dark matter and galaxy 
2- and 3-point correlation functions in the simulation in 
both real and redshift space and extracting the parame- 
ters ci and C2. Since the relation ([9]) holds for the den- 
sity fields on some smoothing scale, we should fit the 
correlation functions for separations comparable to or 
slightly larger than this scale. For the 3PCF, we use 
triangles with u = r 23 /r 12 = 2, r = 9 ft. Mpc, and 
weight equally all configurations with 0° < a < 180° 
using the L200 box. To calculate the likelihood func- 
tion, we use a method s imilar to that described in 
Gaztanag a fe Scoccimarrol (|2005f ). which is based on an 
eigenmode analysis of the covariance matrix. We use 
the jack-knife subsamples in real and redshift space to 
construct covariance matrices, and we use only the dom- 
inant eigemodes with values > y/2/N where N = 16 
is the number of jack-knife subsamples; we found that 
adding further eigenmodes just increases artificially the 
signal. 

The best-fit parameter values from the 3PCF, substi- 
tuted into Eqn. are shown as the short-dashed blue 
curves in Figure [8] Also, in the right panels of Fig- 
ure [2] we show the predicted galaxy 3PCF (solid blue 
curves) using the measured dark matter Qdm and the 
best fitting c\ , c 2 parameters in Eqn. (|10p from the fit at 
r = 9/i _1 Mpc. We see that the agreement with the mea- 
sured galaxy 3PCF is very good in both real and redshift 
space for triangles with r = 9/i _1 Mpc; for configurations 
with r = Gh^ 1 Mpc, the agreement in real space is still 
quite good while in redshift space some deviations in the 
shape-dependence appear. 

In Figure [5] we show the fits for c\ and c 2 for the dif- 
ferent galaxy samples in the L200 box in both real and 
redshift space: all galaxies (top panels), galaxies in the 
absolute range —20 > > —21 (middle), and those in 
the range -21 > > -22 (bottom). The thick oval 
contours indicate the 1 and 2a confidence intervals for a 
Ax 2 distribution with two free parameters, constrained 
using the 3PCF, via Eq. ([1(7)) . and the large points in- 
dicate the maxiumum likelihood values from the 3PCF. 
The thin solid contours show the constraints on c\ and c 2 
using the fit of Eq. ([9]) to the measurements in Figure 
and the vertical line in each panel shows the estimate of 
ci from comparing the 2PCF amplitudes for galaxies and 
dark matter using pairs selected from the triplets used to 
measure the 3PCF. The best-fit values are shown in Ta- 
ble [2] the first two columns show the best-fit parameters 

, using the reduced 3PCF (eq. [TO"]) , the next two 
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Fig. 9. — Constraints on the model bias parameters c\ and C2 in 
real (left panels) and redshift space (right panels), measured in the 
L200 box for all galaxies (top), galaxies with —20 > > —21 
(middle), and galaxies with -21 > > -22 (bottom). Thick 
ellipses correspond to la and 2<r constraints on parameters using 
the reduced 3PCF fit to equation llOl symbols denote the minimum 
X 2 value. Thin ellipses come from fitting the smoothed halo and 
dark matter overdensities (equation [9} shown in Figure [8] The 
vertical lines show the constraint on c\ from comparing the dark 
matter and galaxy 2PCF amplitudes and fitting the ratio with a 
linear bias model, c\ = £, ga i/£, d m- 

columns, cf and c 2 , are calculated using the quadratic 
bias model (eq. [9]) fit to the counts in cells, and in the 
last column c\ is obtained from comparing the 2PCF of 
galaxies and dark matter assuming a linear bias model, 

c l £,gal I dm- 
Ill agreement with previous measurements in surveys 
and simulations, the bias parameters obtained from the 
3PCF are degenerate, resulting in elongated contour el- 
lipses; this could be mitigated to some extent by using 
a larger variety of triangle configurations. In real space 
(left panels of Figure [5]) , we see that the three methods 
of extracting the bias parameters are in rough agree- 
ment: there is a preference for c\ ~ 1 and negative 
c 2 as was found for th e 2dFGRS 3PCF measurements 
(|Gaztanaga et all [20051. Note that the 3PCF fit tends 
to overestimate c\ and to slightly overestimate c 2 com- 
pared to the other two methods. In redshift space (right 
panels), the opposite is true: the 3PCF constraint tends 
to underestimate c\. The right panels of Figs. [5] and 
[9] show that there is a larger discrepancy between the 
3PCF fits and the counts-in-cells fit to Eqn. ([9]) in red- 
shift space, suggesting that the relation pop may not be 
a good representation in redshift space. 

While this comparison with the deterministic local 
bias model is suggestive, to test it more quantitatively 
one should use more triangle configurations and larger- 
volume catalogs to enable more precise calculation of the 
correlation matrices. 

7. SUMMARY 

We have studied the 3PCF of dark matter and galax- 
ies in high-resolution dissipationless cosmological simu- 
lations. The galaxy model associates dark matter ha- 
los with galaxies by matching the halo velocity function, 
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TABLE 2 

Best-fit bias parameters in the L200 box 



Subample 



c Q 



' -2 



1 



All objects r-space 
All objects z-spacc 
—20 < Alp- < ~21 r-space 
—20 < < —21 z-space 
-21 < M^ 1 < -22 r-spacc 
-21 < < -22 z-space 



1.16 
0.86 
1.08 
0.86 
1.42 
0.99 



+0.21 

0.15 
+0.11 
0.14 
+0.28 
0.18 
+0.21 
0.16 
+0.48 
0.29 
+ 0.15 
0-35 



-, + 0.19 

-0.16 

-0 31+° 15 
u - o± -0.11 

-U.^U_ Q 2 3 
-1+0.20 
-0.15 
-+0.33 
-0.35 
-0 20+ ' 21 



-0.19 1 



-0.301 
-0.15" 1 



1.10 
0.97 
1.04 
0.97 
1.21 
1.01 



0.05 
0.12 
+ 0.05 
0.07 
+0.07 
0.16 
+0.06 
0.12 
+0.06 
0.14 
+0.12 
0-05 



-0.20 
-0.32 
-0.08 
-0.34 
-0.36 
-0.43 



+0.19 
0.22 
+0.12 
0.18 
+0.08 
0.19 
+0.22 
0.16 
+0.31 
0.05 
+0.33 
0-07 



1.06±0.09 
1.03±0.08 
1.01±0.09 
1.00±0.08 
1.19±0.10 
1.13±0.13 



including subhalos, with the observed galaxy luminosity 
function by abundance, and has been shown previously to 
provide an excell ent match to observ ed two-point statis- 
tics for galaxies (|Conrov et al.l [20061 . Our primary re- 
sults are as follows: 

1. The reduced real-space 3PCF for both galaxies 
and dark matter has strong dependence on scale 
and shape. The shape dependence of the 3PCF 
strengthens with increasing scale, in agreement 
with previous simulation results for dark matter 
and with expectations from non-linear perturba- 
tion theory. 

2. On small scales, or alternatively with increasing 
time, the shape dependence of Q washes out as 
virial motions within halos replace coherent infall 
on larger scales. 

3. Redshift-space distortions attenuate the shape and 
scale dependence of the reduced 3PCF and weaken 
the evolution with redshift. 

4. The reduced 3PCF shows only weak dependence on 
galaxy luminosity and color; put another way, the 
scaling between the 3PCF amplitude and the 2PCF 
is predicted to be nearly independent of galaxy 
type in this model. The trend of Q with color is 
somewhat stronger than with luminosity: the re- 
duced 3PCF is slightly enhanced for red galaxies 
over blue, especially for elongated configurations. 

5. Our model predictions are in excellent agreement 
with the shape and scale dependence of the galaxy 
3PCF measured in the SDSS when the same bin- 
ning scheme is used. Since the results are highly 
sensitive to the binning scheme, caution must be 
exercised in comparing theory and observations of 
the 3PCF. In combination with earlier results, this 
comparison indicates that a simple scheme in which 
galaxies and dark matter halos and subhalos are 
associated in a one-to-one fashion based on maxi- 
mum circular velocity can provide a good match to 
a wide range of galaxy clustering statistics. 

6. The effect on the 3PCF of changing the selection of 
subhalos (i.e., of connecting galaxy luminosity to 
V ma x,acc instead of V max ,now) is evident on small 
scales (r < 2 /i _1 Mpc) and for elongated configu- 
rations, but is negligible on larger scales. Future 
measurements of the 3PCF will help constrain dif- 
ferent models for the association of galaxy luminos- 
ity and color with subhalo properties. 



7. On scales of order 10/i -1 Mpc, a local, determinis- 
tic bias scheme is in reasonable agreement with the 
galaxy and dark matter distributions of the model. 
The bias parameters extracted from Q ga i are in rea- 
sonable agreement with the 5 ga i-8dm relation in real 
space, less so in redshift space. Nevertheless, the 
redshift-space constraints on the bias parameters 
are in agreement with the 2dFGRS measurements 
of the 3PCF. 
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